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ABSTRACT 

Using a recently proposed nonlinear force-free method designed for single 
vector magnetograms of solar active regions we calculate the instantaneous free 
magnetic energy and relative magnetic helicity budgets in 162 vector magne- 
tograms corresponding to 42 different active regions. We find a statistically 
robust, monotonic correlation between the free magnetic energy and the relative 
magnetic helicity in the studied regions. This correlation implies that magnetic 
helicity, besides free magnetic energy, may be an essential ingredient for major 
solar eruptions. Eruptive active regions appear well segregated from non-eruptive 
ones in both free energy and relative helicity with major (at least M-class) flares 
occurring in active regions with free energy and relative helicity exceeding 4 x 10 31 
erg and 2 x 10 42 Mx 2 , respectively. The helicity threshold agrees well with esti- 
mates of helicity contents of typical coronal mass ejections. 

Subject headings: Sun: chromosphere Sun: corona Sun: flares Sun: surface 
magnetism Sun: photosphere 



1. Introduction 

Active regions (ARs) are formed when considerable, localized magnetic flux emer- 
gence occurs into the solar atmosphere, with characteristic fluxes of the order of 10 22 Mx 



1 Marie Curie Fellow. 



-2 - 



( ISchrijver fc Harveylll994l ). ARs are magnetic structures far from a ground, current-free (po- 
tential) energy state and, as such, they store large amounts of free magnetic energy. Over the 
last two decades, multiple report s highlight the simultaneous accumulatio n of large amounts 



of magnetic helicity in ARs (e.g., lLaBonte et al.ll2007t ISmyrli et al.ll2010l ). Magnetic helicity 



emerges via helical magnetic flux tubes or is being generated by solar differential rotation 
and peculiar photospheric motions and is a way to quantify the stress and distortion of the 
magnetic field compared to its potential-energy state. The free magnetic energy release is 
fragmented in ARs and fuels solar flares and/or coronal mass ejections (CMEs) that tend to 
relax the magnetic configuration. Co ntrary to mag netic energy, helicity cannot be efficiently 
removed by magnetic reconnection (jBergerl Il984l ) . If it is not transferred to larger scales 
i n the Sun yia existing m agnetic connections, it can only be expelled in the form of CMEs 
( ILowl Il994j ; iDeVord |2000| ). As a result, an isolated, confined magnetic configuration with 
accumulated magnetic helicity cannot relax to a potential-field configuration. 

Although the importance of storag e and release of free magnetic energy in ARs is widely 
acknowledged for solar eruptions (e.g., ISchrijverl 120091 ). the role of magnetic helicity is still 
under de bate. Indeed, it has been demonstrated tha t helicity is not necessary for eruptions 
to occur (IPhillips et al.ll2005l ; IZuccarello et al.l 120091 ) . On the other hand, observational and 
modeling works have shown that ARs with large dominant left- or right-handed helicity give 



rise to more and/or major erupt i ons fe.g.. iNindos fe Andrews! 120041 ; iTorok fc Kliem 



2005 



LaBonte et al.ll2007t lNindosll2009l ; iGeorgoulis et al.ll2009h . Instrumental to this debate is the 
lack of robust methods to calculate the (relative to a reference potential field) instantaneous 
magnetic helicity budget in ARs. Existing methods basically re strict to either integrating 
in time the relati ye helicity injection rat e (IBerger fc Field! 119841 ) or evaluating the relative 
helicity formula ( jFinn fc Antonsenl Il985l ; iBergerl Il999l ) in a volume by means of a three- 
dimensional magnetic field extrapolation and a respective gauge- dependent expression for 
the vector potential. The helicity injection rate depen ds on the photosphe ric velocity field 
whose inference involves significant uncertainties (e.g.. IWelsch et al.l 120071 ) while the "vol- 
ume" helicity calculation depends on the model-depende nt nonlinear force-free field extrapo- 
lation , also subject to uncertainties and ambiguities (e.g., 
20081 . and references therein). 



Schrijver et al.ll2006l ; iMetcalf et al. 



Recently, IGeorgoulis et al.l ((2012J) (hereafter GTR12) proposed a general force-free ap- 
proach to self-consistently calculate the instantaneous magnetic free energy and relative 
helicity budgets from (photospheric or chromospheric) vector magnetograms of solar ARs. 
The method does not rely on any magnetic field extrapolation but it uses a magnetic connec- 
tivity matrix that may be inferred by extrapolations, among other methods. As in GTR12, 
we use in this work a simulated annealing method that converges to a unique connectivity 
matrix to calculate the free magnetic energy and relative magnetic helicity for a large sam- 
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pie of active-region vector magnetograms, seeking a statistically robust correlation between 
these physical parameters, if any. Section 2 briefly describes the method, its application and 
results are discussed in Section 3, while Section 4 summarizes our findings. 



2. Nonlinear force-free magnetic energy and helicity budgets in solar ARs 



In GTR12 we combined the linear force-free technique of iGeorgoulis fc LaBontd (120071 ) 
with the properties of the mutual helicity as discussed by iDemoulin et al. to derive 

the instantaneous nonlinear force-free (NLFF) field energy and helicity budgets in an AR. 
In particular, a vector magnetogram of a given AR is translated into a collection of slender 
force-free flux tubes with known footpoints, flux contents, and variable force-free parame- 
ters. These flux tubes are defined after determining the magnetic connectivity matrix in 
the AR's magnetogram, which provides the flux contents committed to the connection be- 
tween opposit e-polarity flux parti ti ons. The method of choice to calculate this matrix was 
introduced by IGeorgoulis fc Rustl (120071 ) and was revised by GTR12. It uses a simulated 
annealing method that globally (within the field of view) minimizes the connection lengths, 
at the same time guaranteeing that only opposite-polarity flux elements will be connected to 
each other. This approach emphasizes complex active regions with intense magnetic polarity 
inversion lines (PILs) in their lower boundaries. As explained in GTR12, the method gives 
unique results, contrary to model-dependent NLFF field extrapolations. 

For a collection of N slender flux tubes the free magnetic energy E c is the sum of a self 
term E r ... due to the internal twist and writhe of each flux tube, and a mutual term E c ,, 

'-self " " t-mut " 

due to interactions between different flux tubes. GTR12 provide a lower limit of the free 
energy E c for a given connectivity that (i) assumes no winding of a given flux tube around 
others, and (ii) neglects contributions from potential flux tubes induced by the existing tubes 
in a space-filling, force-free magnetic configuration. This expression for E c reads 



E c setf + E, 



C-mut 



N N N 

Z=l 1=1 m=l,l^m 



(1) 



In Equation ([!]) A and S are known fitting constants, d is the pixel size of the magnetogram 
and and a\ are the respective flux and force- free parameter of flux tube I. Cf™ h is the 
mutual-helicity factor of two arch-li ke (not winding around each other) flux tubes. Inference 
of this factor was first discussed by lDemoulin et al.l (120061 ) and was later refined by GTR12. 
Refinement included the case of a "matching" photospheric footpoint in a pair of flux tubes. 
"Matching" footpoint means that the like-polarity footpoints of a given flux-tube pair are 
within the same magnetic partition; therefore, they are unresolved (considered as coinciding) 
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by the method. To calculate Cf™ one nee ds the relative footpoint locations of flux tubes / 
and m and an assessment of whether I is "above" m, or vice versa. For each of these cases one 
calculates a different £^ c?l -value in cases of "matching" footpoint and intersecting footpoint 
segments for the flux-tube pair. By definition, these two values have opposite signs. In case 
of non-intersecting segments, the two possible £^ c/l -values collapse to a single value. Cf^ h is 
always a real number with an absolute value lower than one. The selected Cf^ h value is the 
one that assigns a positive increment to the overall free magnetic energy because of the flux- 
tube pair. For a pair (l,m) of flux tubes with force-free parameters on and a m , respectively, 
this means (a; + a m )£f™ h > 0. 0\m h is assumed equal to zero if the respective energy 
increment can only be negative (this applies exclusively to the non-intersecting-segments 
case), which is not a physical solution. 

The respective relative magnetic helicity H m for the collection of N slender flux tubes 
is also the sum of a self (H mself ) and a mutual (H mmut ) term. GTR12 derive 

N N N 

H m = H mselI + H mmut = 8nd 2 A ^ a&f + ^ E • (2) 

1=1 1=1 m=l,lj^m 

A detailed derivation of uncertainties for both E c and H m is also provided in GTR12. 

3. Application to observed solar active region magnetic fields 

3.1. Data selection 



We have collected an extended sample of 162 photospheric and low chromospheric 
vector magnetograms obtained by the g round-based Imaging Vector Magnetograph (IVM; 



Mickey et al.lll996t lLaBonte et al.lll999l ) of the University of Hawaii's Me es Solar Observa- 



tory and by the space-based Spectropolarimeter (SP; see description in iLites et al.l 120081 ) 
of the Solar Optical Telescope (SOT) onboard Hinode. The IVM provides complete Stokes 
profiles of the Fe I 630.2 5 nm photospheric lin e (for earlier data) and of the Na I D 1 589.60 
nm chromospheric line (ILeka fe MetcalJ 120031 1 . for data recorded after 2003, with a spatial 
sampling of 0.55 arcsec per pixel (full) or 1.1 arcsec per pixel (binned). The SP provides 
full Stokes profiles of the Fe I 630.25/630.15 nm lines with a maximum spatial sampling of 
0.16 arcsec per pixel. IVM magnetogram edges were carefully cropped to remove instru- 
mental border artifacts. To remove the intrinsic 180° azimuthal ambiguity in the vector 
magne tograms we a pplied the non-potent ial field calculation (NPFC) method of iGeorgoulis 
(120051 ). as revised in iMetcalf et al.l (120061 ). As typical uncertainties for the line-of-sight and 
transverse field components we used (SBi, 5B tr ) = (50, 100) G, for IVM data, and (8B h 5B tr 
) = (5, 50) G, for SOT/SP data. 
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Our sample of 162 vector magnetograms corresponds to 42 different ARs, distinguished 
between "non-flaring" and "flaring" ones. The first category involves 18 ARs that have 
not hosted a higher than C-class flare, while the second involves 24 ARs with at least one 
M-class flare. Non-flaring ARs include NOAA ARs 8844, 9114, 9635, 9845, 10050, 10254, 
10323, 10349, 10536, 10939, 10940, 10953, 10955, 10956, 10961, 10963, 10971, and 10978. 
Flaring ARs include NOAA ARs 8210, 9026, 9165, 9393, 9415, 9632, 9661, 9684, 9704, 9773, 
10030, 10039, 10162, 10365, 10375, 10386, 10484, 10488, 10501, 10570, 10596, 10656, 10930, 
and 10960. For a few ARs we have timeseries of vector magnetograms spanning from a few 
hours of IVM data (e.g., NOAA ARs 8844, 9165) to a few days of SOT/SP data (e.g., NOAA 
ARs 10930, 10956). Moreover, our sample includes both short-lived emerging flux regions 
(i.e., NOAA AR 8844) and persistent, large and complex regions (e.g., NOAA ARs 10488, 
10930). 

Figure [1] shows the central heliographic position of the selected 162 magnetograms. Lo- 
cations cover a latitudinal area of ±25° and a meridional zone of —40° to +60° . Larger 
central-meridian distances were avoided to avoid extreme projection effects. The calcula- 
tions below involve the local, heliographic field components on the image plane of each 
magnetogram. 

3.2. Magnetic energy — relative magnetic helicity budgets and their relation 

Using equations (p]) and (|2J) we calculate the free magnetic energy and relative magnetic 
helicity budgets for the selected 162 vector magnetograms. The results suffice to construct 
the free-energy - relative helicity diagram (hereafter energy - helicity [EH] diagram) of solar 
ARs, shown in Figure [2J Two main conclusions stem from this plot: first, there is a nearly 
monotonic dependence of direct proportionality between E c and H m ; flaring ARs tend to 
show both large free energies and large amplitudes of relative helicity. Second, flaring and 
non-flaring ARs appear well segregated; excesses of ~ 4 x 10 31 erg in free magnetic energy 
and of ~ 2 x 10 42 Mx 2 in relative magnetic helicity tend to bring ARs into major flaring (or 
eruptive, in general), territory. Within this area, there is no appreciable segregation between 
M- and X-flaring ARs. Later in this Section we discuss the (few) exceptions to the above 
general assessment. A similar diagram (not shown) between the free magnetic energy and 
the unsigned (total) magnetic flux shows significantly less segregation between flaring and 
non-flaring ARs than that of the EH diagram. 

The least-squares best fit between \H m \ and E c reveals a scaling of the form 

97 4^ 

log \H m \ oc 53.4 - 0.0524 (log£ c ) a653 exp — — , (3) 

log E c 
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while, in general, a simpler logarithmic scaling of the form 

\H m \ oc 1.37 x 10 14 £ c a897 (4) 

also works quite well. Both have a goodness of fit ~ 0.7 with the Kolmogorov-Smirnov test 
giving also similar significance levels. Equation (|3]) implies that relative magnetic helicity 
builds up with faster fractional rates than magnetic free energy in the initial stages of active 
regions, however this is not conclusively supported by the derived fit. 

Figure [3] shows the calculated uncertainties for the helicity magnitude of the studied 
162 magnetograms, ordered by increasing free magnetic energy. The respective uncertainties 
for the free magnetic energy are not shown because they are quite smaller: below the free 
energy threshold the mean error is ~22% while above it the mean error becomes ~7%; hence, 
they do not affect significantly the EH diagram of Figured especially its high-energy part. 
Uncertainties in the relative helicity are higher: below the relative helicity threshold the 
mean error is ~50% while above it the mean error drops to ~14%. As with free-energy 
uncertainties, the respective helicity uncertainties do not alter significantly the EH diagram 
of Figure [2J Point taken, there are some notable exceptions (see below). 

Before discussing exceptions, let us briefly discuss how the two thresholds in free energy 
and relative helicity divide the EH diagram of Figure [2j Regions a and c include populations 
of small E c / small \H m \ and large E c / large \H m \, respectively and hold the vast majority of 
points in the diagram, thus reflecting the nearly monotonic free-energy - helicity dependence. 
Region b indicates ARs with a large free energy but with a small, or relatively small, relative 
helicity budget. Given that magnetic helicity is a signed quantity, a population in Region 
b might also include ARs with significant, but similar, amounts of both senses of helicity. 
Other than four (4) magnetograms of X-flaring ARs and one (1) magnetogram of a M-flaring 
AR (discussed below) there is no clear population in Region b. This is important evidence 
to the existence of a dominant sense of helicity in strongly helical ARs. Finally, Region 
d indicates ARs with small, or relatively small, free energy and a large relative helicity 
budget. A population in this Region would be troublesome to interpret or would indicate 
a problematic numerical approach to calculate E c and \H m \. Fortunately, other than some 
typical scatter in the EH diagram, close to both thresholds, there is no population in Region 
d. 

We now discuss some notable exceptions to the physical interpretation of the EH diagram 
of Figure [2J first, four magnetograms of flaring NOAA AR 9165 in Region d are within 
uncertainties from the energy threshold, hence, they are not considered exceptions. In Region 
c there are 7 magnetograms of non-flaring ARs (NOAA ARs 10963, 9845, 10953, 10978, 
10536, 10323 and 10349, by increasing free magnetic energy) that clearly go well into major- 
flaring territory. Most of these ARs share a common feature: rather than hosting one 
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or more major flares (they were apparently capable for this based on their free magnetic 
energy and helicity budgets) they released substantial free energy in lengthy series of C-class 
flares. For example, NOAA AR 10536 produced 34 C-class flares within two weeks (the 
highest being a C7.7 flare) while NOAA AR 10349 produced 36 C-class flares in ten days. 
Two magnetograms (NOAA ARs 10978 and 10323) where recorded after a C-class flare. In 
Region b there are 4 magnetograms of X-flaring ARs (NOAA ARs 10030 and 10930) and 
one magnetogram of a M-flaring AR (NOAA AR 10960) that show deficits in their relative 
helicity budgets. All X-flaring ARs (marked in Figure [3]) show relatively high uncertainties 
in relative helicity, with fractional uncertainties in the range of 0.5 to 4.7, capable of placing 
them in Region c. Moreover, and related to the large uncertainties, we cannot rule out 
possible white-light contamination due to flaring in these regions: data for the magnetogram 
of NOAA AR 10030 were recorded right after two X and M-class flares while magnetograms 
of NOAA AR 10930 were recorded between swarms of B- and C-class flares. For the M- 
flaring NOAA AR 10960 we cannot exclude the possibility of similar helicity budgets of both 
senses. We cannot assume the same for previously discussed NOAA ARs 10030 and 10930 
because they participate in the sample with timeseries of magnetograms that, excluding the 
ones in Region b, reside exclusively in Region c. Finally there is one magnetogram of a 
M-flaring (NOAA AR 10501) and two of X-flaring (NOAA ARs 10375 and 10386) ARs that 
fall short of the energy and helicity thresholds, residing into Region a. Significant errors 
(both in free magnetic energy and relative helicity) or contamination by white-light flare 
emission may also be responsible for this effect (the NOAA AR 10375 magnetogram was 
recorded during the declining phase of a C2.5 flare) but, moreover, it is kno wn that intense. 



local i zed magnetic flux eme rgence can also force major eruptive flares (e.g.. lNitta fc Hudson 



200ll ; IZhang fc Wangj|2002l ) before an AR manages to build a strong PIL and, consequently, 



large budgets of E c and H m . 

To test the validity of our free magnetic energies E c , we compare them with the lowest 
possible free magnetic energy E CWT that corresponds to a given amount of relative helicity 
for the NLFF field. This limiting value is th e LFF free magnetic energy corr esponding to 



this helicity, per the Woltjer-Taylor theorem ( IWoltierlll958l : lTavlorlll974 Il986f ). As GTR12 

have shown, E CWT is given by 

Hi 
(87rd)M$ 2 

All calculated free magnetic energies E c must be larger than this value. As Figure H] clearly 
demonstrates, this condition, within applicable errors, is always satisfied. 



77 __ 11 m /c-\ 

CWT fa^r1\2 Aif>2\ ■ 
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Conclusions 



We have applied a new NLFF field method to calculate the instantaneous free magnetic 
energy and relative magnetic helicity budgets of solar active regions. On a sample of 162 
such magnetograms we report, for the first time, (1) a nearly monotonic relation of direct 
proportionality between the free magnetic energy and the relative magnetic helicity in ARs 
and (2) the existence of thresholds ~ 4 x 10 31 erg and ~ 2 x 10 42 Mx 2 , for the free energy 
and the relative helicity, respectively, for ARs to host major, typically eruptive, flares. As 
magnetic helicity is a signed quantity, the monotonic dependence between it and the free 
energy implies that, in spite of simulations reporting that helicity is not necessary for so- 
lar eruptions (see Introduction), real flaring/eruptive ARs show a significant budget and a 
dominant sense of magnetic helicity. This finding appears to suggest an important role of 
magnetic helicity in solar eruptions. However, the details of this role are yet to be uncovered. 

An additional important finding refers to the inferred threshold values for E c and \H m \. 
While a thresh old of 4 x 10 31 erg in free energy suffices to justify at least one M-class flare 
( ]Hudsonll201ll ). a threshold of 2 x 10 42 Mx 2 in relative helicity is in excellent agreement with 
estimat ed helicity budgets of typ ical CMEs (~ 2x 10 42 Mx 2 bv lDeVore! feoooh : (1.8-7) x 10 42 
Mx 2 by iGeorgoulis et al.l (120091 )). This may indicate that active-region CMEs occur only 
when source ARs can shed the required helicity budget for them, in line with theoretical 
CME interpretations as means of relieving the Sun from its excess helicity. 

The above results are a first step to assess and uncover the possible combined role of 
free magnetic energy and relative magnetic helicity in solar eruptions. We do not envision 
or expect use of these results for eruption forecasting purposes at this stage. However, these 
findings instigate focused research on the detailed, quantitative role of magnetic helicity 
in solar eruptions that could even lead to the formulation of credible eruption initiation 
mechanisms. We intend to investigate such possible physical links in the future. 
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Fig. 1. — Central heliographic positions for our 162 vector magnetograms. Blue diamonds 
correspond to non-flaring (up to C-class flaring) ARs while red squares and asterisks corre- 
spond to M- and X-class flaring ARs, respectively. 
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Fig. 2. — The free energy - relative helicity diagram of solar ARs. Blue diamonds, red squares 
and red asterisks correspond to non-flaring, M- and X-class flaring ARs, respectively. Dashed 
lines indicate the estimated thresholds for relative magnetic helicity (~ 2 x 10 42 Mx 2 ) and 
free magnetic energy (~ 4 x 10 31 erg) above which ARs seem to give major flares almost 
exclusively. These thresholds divide the diagram in four Regions, labeled a, b, c and d (see 
text). The dotted and dash-dotted lines denote the least-squares best fit (Equation ([3])) and 
the least-squares best logarithmic fit (Equation fll])), respectively. Some NOAA AR numbers 
are also indicated in Figure Eland are discussed in the text. 
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Fig. 3. — Uncertainties for the difference between the relative helicity magnitudes and the 
inferred relative helicity threshold of 2 x 10 42 Mx 2 , ordered by increasing free magnetic 
energy The dashed vertical line separates ARs below (left) and above (right) the respective 
free magnetic energy threshold. Symbols and marked NOAA AR numbers are the same with 
those of Figure EJ 
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Fig. 4. — Scatter plot for free magnetic energies E c and their respective Woltjer- Taylor 
minima E CwT (Equation (jSJ)). The dotted line denotes equality between the two energy 
estimates. 



